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Abstract 

Diffraction orders in the continuous wave regime generated by a Ronchi transmission grating 






C^ ■ in a standard thresliold configuration are sliown theoretically to violate quantum duality for a 

O^i locally real representation. The phenomenon superficially resembles Rayleigh anomalies but is 

notably distinguished from those anomalies by a prediction of probability non-conservation. This 

\jQ ' prediction is experimentally tested with a 633 nm laser beam at normal incidence on gratings giving 

o: . . 

en ■ that threshold condition for the ±3*" order pair. Transient intersection of the order with an 



o 



% 



independent 633 nm laser beam demonstrates a duality-violating probability non-conservation in 
good agreement with the theoretical prediction. 



I. INTRODUCTION 

The compact mathematical formahsm of quantum mechanics developed in the late 1920's 
continues as an extraordinarily successful representation of physical phenomena. The very 
compactness itself is widely perceived as evidence of the validity, beauty, and completeness 
of that formalism. This perception has persisted despite the evolving understanding over 
the years that strict adherence to the formalism imposes fundamental violations of classical 
physical reality. The resolution of these violations, given the constraints of that strict 
adherence, was understood by Bohr and others to necessarily impose a non-real and non- 
local representation of the physical world. This representation is often referred to as the 
probabilistic interpretation of quantum mechanics (PIQM) and consists of those postulates 
thought to rationally express the physical implications of the underlying compact quantum 
formalism. 

The predictive validity of the quantum formalism is widely acknowledged. Nevertheless, 
the PIQM departures from local realism are not universally accepted. Lepore and Selleri 
[1] contend that "The development of local realism since the 1935 Einstein, Podolsky, and 



Rosen's paper [2| is by far the most profound criticism of quantum mechanics..." in reference 
to the probabilistic interpretation. Popper raises compelling philosophical arguments in 
questioning the validity of P1QM.|3|] 

The notable non-local properties of PIQM are manifestations of entanglement and wave- 
particle duality. Any viable local alternative to PIQM must minimally demonstrate an 
alternative self-consistent basis for both of these phenomena. In this regard, one of the 
present authors proposed a locally real representation of quantum mechanics that gives 
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agreement with performed experiments for correlated photons and particles j^ and is not 
restricted by Bell's theorem. |5| The other phenomenon, wave-particle duality, is examined 
by the authors in ref. |6|. In this report we again examine wave-particle duality and propose 
a locally real representation of quantum mechanics, identified here as LRQM, that retains 
the accepted underlying quantum formalism with minimal modification. 

As a preliminary matter, we first consider some elementary aspects of photon phenomena 
from the particular perspective of LRQM as distinguished from PIQM. In the discrete regime, 
a single photon is a physically real wave packet structure on which a real energy quantum 
resides. The LRQM wave function $ characterizes the wave amplitude of this structure but 



is not itself a physical representation of any resident energy quantum. The separability of 
;hese entities is identified with de Broglie's initial representations of quantum phenomena 
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and is intrinsic to many locally real representations. 

Born interpreted the squared modulus of the wave function evaluated at particular co- 
ordinates as a measure of relative probability of finding a particle on a wave packet. [8|] 
This initial restricted version of "Born's rule", applied here to electromagnetic radiation for 
LRQM, is critical in its treatment of the squared modulus |$| as a relative and not an 
absolute probability flux density. Saxon reminds us of the fundamental utility of the un- 
normalized I'I'I as the sufficient determining factor of relative positional probability. J9| The 
subsequent, commonly accepted interpretation of Born's rule incorporates normalization of 
the wave function giving unit absolute probability for the squared modulus integrated over 
the entire wave packet. Normalization is a seemingly natural and rational modification of the 
restricted version and was fully consistent with the then developing formulation of PIQM. 
This modification provides PIQM with a linkage of particle-like and wave-like properties 
that compactly expresses both in a wave function $P7 but necessitates the characteristic 
interpretations of duality and entanglement for particular phenomena. 

Certainly, for an ordinary photon the resident particle-like energy quantum has unit 
absolute probability of existing somewhere on the wave packet and normalization of the 
LRQM wave function $ in this case serves as a mathematical convenience in equating a 
unit-valued expression of energy conservation and a unit- valued probability represented by 
the integrated square modulus |$| . Nevertheless, in LRQM an initial normalization of 
$ must not be arbitrarily re-applied to an evolving wave function since processes such 
as interference occurring during that evolution may create or annihilate wave amplitude. 
Effectively, LRQM separates a purely wave-like $ from a PIQM $p/ and provides that $ 
with the degree of freedom to scale independent of resident quanta. We continue to use the 
term "probability" here in LRQM bearing in mind that its usage is potentially misleading 
since that term suggests equivalence to a mathematical absolute probability, an equivalence 
manifested in PIQM as duality. 

Classically, the probability flux density |$| is recognized as a wave intensity. For a 
discrete photon incident on an idealized beam splitter, the wave intensity fractionally divides 
onto the output channels in accordance with the transmission and reflection coefficients of 
the beam splitter. In the present example, we choose for convenience a 50:50 beam splitter 



for which the output packets are similar to the incident packet but with half the intensity 
giving each a relative probability of P = 0.5 when the incident packet is assigned P = 1. 
Then, when the energy quantum on that incident packet reaches the beam splitter, the 
quantum randomly transfers onto either of the emerging outgoing P = 0.5 packets with 
equal probability. In the discrete regime this implies that for each incident photon, one 
of the outputs is an "empty" wave that is totally "depleted" in energy quanta relative to 
its wave packet probability P = 0.5. Conversely, the other output is "enriched" in energy 
quanta relative to probability in the regard that the single energy quantum resides on a wave 
packet that is now P = 0.5 instead of P = 1. The prediction of empty waves in the discrete 
regime has been the subject of many investigations seeking to differentially test local realism 



and PIQM as in a series of papers by Croca et al. lu\ as well as in numerous others such as 
refs. pjj]. These investigations, frequently using a beam splitter to generate an empty wave, 
are necessarily restricted to the very weak wave intensities associated with discrete photon 
beams. 

As we move from the discrete photon regime to a multiple photon beam in the continuous 
wave (cw) regime, the relevant total wave function is conventionally constructed from a 
summation over amplitudes of the constituent wave packets. (Our particular interest here 
is mono-energetic coherent beams.) The wave function, which we continue to identify as $, 
is used to express the beam's intensity |$| , again a probability flux density. Integration of 
I'I'I over some selected beam segment gives an inclusive probability P. The corresponding 
inclusive quanta in that segment have a total energy E. We can set P = E = 1 in arbitrary 
units for the beam segment. 

When this beam is incident on a 50:50 beam splitter, we again have P = 0.5 on each of the 
output channels. However, random statistical distribution of the inclusive quanta onto the 
output channels yields E vanishingly close to 0.5 on both channels as the inclusive number 
of quanta becomes large. An immediate consequence of the cw regime for the beam splitter 
is that testability for empty waves is no longer feasible. Then, although local realism does 
not expressly prohibit empty waves in the cw regime, the inherent statistical distribution 
process is intuitively expected to restrict all mechanisms for generating empty waves to the 
discrete regime where conclusive experimental verification is marginalized. 

A mechanism that alters the proportionality of probability and energy on a photon beam 
in the cw regime is seemingly as improbable as Maxwell's hypothesized mechanism for 



selectively sorting particles based upon their respective kinetic energies. [12 1 

In Section III. we present the theoretical basis for a mechanism that, from the perspective 
of LRQM, is predicted to selectively alter the proportionality of probability and energy. In 
Section IV. we report on the experimental realization of this phenomenon. 

II. BACKGROUND 

The theoretical basis developed here for the LRQM wave-like component is substantially 
classical. This basis does not obviate the underlying quantum mechanical formalism such 
as the wave function associated with quantization of the electromagnetic field, but it does 
modify the scaling of a separable, purely wave-like component $. As a consequence, the 
PIQM duality proportionality of wave probability (based on that $) and particle-like en- 
ergy quantum is potentially violated. A mechanism for this duality violation, trivially but 
transiently realized in the discrete regime for LRQM, is however not immediately obvious in 
the cw regime. Before beginning the examination of mechanisms that achieve this "duality 
modulation" in the cw regime, we return to the example of a beam splitter once again in 
order to develop some basic operational definitions and principles relevant to LRQM. 

When a photon is incident on a beam splitter, the transfer process of the energy quantum 
to one of the relative probability wave packets on the output channels is itself of fundamental 
significance to LRQM. For illustrative purposes we again choose a 50:50 beam splitter. The 
quantum on an incident P = 1 wave packet enters a zone at the face of the beam splitter 
where the quantum randomly transfers onto one of the two emergent P = 0.5 output wave 
packets. Unlike their treatment in PIQM, those output wave packets remain as real entities 
at that probability value on both output channels and one packet does not undergo a collapse 
when a measurement is made on the other packet. 

From the perspective of LRQM, the transfer at the beam splitter must be treated purely 
as that of the energy quantum. Wave packet probability distribution is a non-quantized 
deterministic process at a device such as the beam splitter. The energy quantum arriving at 
the beam splitter surface randomly transfers onto one of the two emergent outputs causally 
mediated by their relative probability distribution consistent with Born's rule. |8| The wave 
structures of the two output packets are otherwise unaltered by the random presence or 
absence of the quantum on a particular output packet. (Conversely, we will subsequently 



consider full photon transfers in other contexts that are consistent with both PIQM and 
LRQM, i.e. duality is not violated. Such transfers may involve scattering conditions that 
alter the trajectories of incident wave packets which are then accompanied by the propor- 
tionate energy quanta that had resided on the incident packets.) 

In the interests of formalizing the proportionality between wave packet probability and 
energy quantum, we assign a value to their proportionality. Each photon on a discrete 
beam is represented by an "occupation" value Q = E/P defined as the ratio of the energy 
quantum and the total probability of its wave packet, [g] As a baseline reference, we consider 
"ordinary" photons that might be generated by common atomic emission processes. We 
assign E = 1 io the energy quantum present on a wave packet of probability P = 1 giving 
fi = 1/1 = 1 for these ordinary photons in arbitrary dimensionless units. 

If a discrete beam composed of such ordinary photons is incident on the 50:50 beam 
splitter, the output packets have Vt = 0/0.5 = and Q = 1/0.5 = 2. Any fi < 1 signifies 
that the outgoing wave packet is reduced in energy relative to its associated wave packet 
probability and is referred to as "depleted" . In the extreme case oiQ = 0, the wave packet is 



totally depleted and is appropriately referred to as an empty wave.|lll| For i7 > 1, the wave 
packet is said to be "enriched" . In the present case oi Q = 2, the single quantum resides on 
a P = 0.5 wave packet compared to the P = 1 wave packet of the incident photon. 

In principle, interference of either outgoing wave packet P = 0.5 with an independent 
photon beam would result in the same visibility, unaffected by the presence or absence of 
the quantum on that outgoing wave packet. Fundamentally, a measurement such as interfer- 
ence visibility assesses a beam's wave-like property whereas a direct detector measurement 
assesses a beam's particle-like (energy quantum) property. 

As a matter of practical consideration, measurements of the wave-like property in the 
discrete regime, as provided for example by beam splitter outputs, are experimentally prob- 



lematic. [10| Accordingly, we are motivated to seek a duality violation in the cw regime 
where experimental verification is significantly enhanced. 

In the cw regime, the irradiance I is used to describe the particle-like component of the 
photon beam where, conventionally, the formal units are those of an energy flux density. 
In contrast, the corresponding wave-like component is identified as the intensity VT = |$| 
which, in the LRQM context, is exclusive of the energy content of the wave. In analogy to 
our use of arbitrary units for the discrete beam, we are free to set unit values for I and W 



on an ordinary incident cw beam. Any value specification of flux densities such as I and W 
implicitly refers to a sample point on the beam and typically that sample point gives the 
maxima of those values e.g. at the beam centroid for a Gaussian cross-section. 

With a specified / and W on an incident beam we have Vt = I/W. Equilibration of 
energy quanta ensures that the Q proportionality is maintained throughout the beam. Then 

w 
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where the areal integration is over the beam's cross-section and the integrands are the 
respective coordinate-dependent values of / and W . An inclusive energy E and probability 
P on a selected beam segment may be selected by a choice of temporal integration limits 
spanning some selected At. For that At, we are free to assign arbitrary units to these 
quantities on the incident beam that give E = P = 1 whereby Q = 1. Clearly, from the 
perspective of LRQM, fi = 1 is maintained on the outputs where, for a 50:50 beam splitter, 
P = 0.5 and, because of random statistical equilibration, E = 0.5. 

As a result of the common equivalence of the I/W and E/P ratios on any given beam 
in the cw regime, we will have frequent occasion to consider either the /, W (maximum flux 
density) quantities or the E, P (integrated) quantities as a matter of illustrative convenience. 
Within the context of the present analysis, when E, P is the more convenient pair, the energy 
flux (power) AE/At and the probability flux AP/At would be equally convenient. However, 
we choose E, P in the interests of continuity with our earlier discussion of discrete photons. 

As a reference baseline, an incident beam consisting of photons generated from conven- 
tional sources is identified as ordinary with Q = Ihj appropriate assignment of the arbitrary 
energy and probability units. Then the identifications of a cw beam as depleted or enriched 
are equivalently expressed by a respective Q < 1 ot Q > 1. 

We consider one last example involving the beam splitter, again 50:50, in which both 
input channels have mutually coherent similar incident cw beams with intensities 

Wih = Wi, = 0.5 (2) 

directed at the same point on the beam splitter and aligned so as to generate spatially 
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FIG. 1: Amplitude moduli associated with two identical mutually coherent beams coincident on 
a beam splitter depicted with lateral displacement for illustrative clarity. Outputs, that form 
resultants, continue propagating as real entities. For destructive phasing of Aqw and Aohv^ the 
two incident beams are fully expressed by the resultant A^h as an intensity output W^h- Arbitrary 
numerical values are assigned for comparative purposes. 

coincident outputs as identified in Fig. 1. In terms of the respective probabilities, 



ih 



P,., = 0.5. 



(3) 



In this current example, we continue the convention of designating beam quantities as 
"incident" prior to interaction with a mechanism such as the beam splitter and "output" 
for the prompt post-interaction manifestations of those quantities. These post-interaction 
output quantities then generate physically distinguishable quantities by interference that 
are designated as "resultant" . Clearly this distinction of output and resultant designations 
was superfluous for the previous example of a single beam incident on a beam splitter. 
Nevertheless, the utility of these designation distinctions will be evident when we consider 
duality modulating mechanisms. 

In the present case, the two incident beams intersect at a common point on a beam 
splitter. In Fig. 1, the beams are laterally displaced from that common point in order 

8 



to clearly depict the various output and resultant components emanating from that point. 
The amplitude moduli of the two reflected and the two transmitted output beams all have 
identical values 

^ohv = ^ovh = ^ohh = ^ovv = 0.5 (4 J 

with corresponding common intensity values 

Wohv = Wovh = Wohh = Wo.. = 0.25 (5) 

or, equivalently in terms of respective probabilities, 

Then for total probabilities in the transition from incidence to output, 

P^ = Po = 1, (7) 

and probability is conserved in this incident— T-output transition. 

Interference between the output beam amplitudes generates a pair of resultant beams 
with moduli dependent upon relative incident beam phasing. For our purposes here, we 
choose a particular phasing that gives amplitude moduh 

Arh = 1 (8) 

and 

Ar, = (9) 

with wave intensities 

Wrh = 1 (10) 

and 

Wr, = 0. (11) 

The output pair Aohh and Aovh is fully contiguous to its associated resultant Arh, sharing 
a common origin at the beam splitter surface. The same applies to Aom, A^hv and A.f.v 
In both cases, the members of either output pair are not separately experimentally mea- 
surable because interference of those members at their common origin promptly yields the 
experimentally measurable resultant wave. Nevertheless, these considerations do not alter 



the continued reality of the output wave structures as they propagate contiguous to and are 
physically expressed by their respective resultant waves. 
Integration yields the associated probabilities 

Prh= fWrh{r,t)dadt = l (12) 

and 

P^„ = / Wrv{r,t)dadt = (13) 

where again, by appropriate choice of integration limits in arbitrary units, each has the 
same numerical value as that of the corresponding (maximum) wave intensity. With the 
total resultant probability 

we have conservation of probability in the output— T-resultant transition as well as in the 
incident— )-output transition. 

For our particular choice of relative phase, Prv is identically zero and Prh = 1- (That 
choice could also have been generalized as any relative phase but in anticipation of a rele- 
vant analogy in the next section, we choose a phase that gives a null value for one of the 
resultants.) For Prv, the constituent output waves continue to propagate as real but n out- 
of-phase entities with a null sum. Similarly, the constituent in-phase output waves generate 
a resultant Prh = 1 that exceeds the sum of those output probabilities. 



ohh 



+ Po,h = 0.5. (15) 



These observations are very elementary and would have been superfluous for the absolute 
mathematical probabilities of PIQM. However, in the context of the real wave structures of 
LRQM, there is a significant process that should be emphasized here. Interference, which 
nullifies the resultant of the two output probabilities on one channel, fortuitously "amplifies" 
the output sum probability on the other channel to precisely compensate for that nullification 
and conserve probability in the output— )-resultant transition. The criticality of interference 
in conserving probability on these real wave structures suggests that a disruption of this 
precise compensation by interference may provide a mechanism for duality modulation. 

Then with fi's computed from either the /, W or the E, P pairs, 
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respectively for the horizontal, vertical, and total resultant beam occupation values. These 
LRQM values are fully consistent with PIQM duality. 

The various beam splitter phenomena treated here from the perspective of LRQM are, 
nevertheless, almost universally represented in the literature by a formalism consistent with 
PIQM. However, in the next section we analogously apply LRQM to particular grating sys- 
tems and demonstrate the basis for a duality-violating mechanism that is not representable 
by PIQM. 

III. THEORETICAL BASIS FOR DUALITY VIOLATION 

A. Gratings with dense sampling 

The general mechanism we consider here for modulating the ratio Q of beam energy 
and probability is a conventional transmissive "picket-fence" grating consisting of a regular 
linear array of parallel opaque bands with a periodicity p forming intervening slits of width 
w. Typically for optical gratings of this type, the opaque bands are thin metallic depositions 
on one side of a transmissive substrate. In our analysis, we treat that side as the exit face of 
the transmission grating with an incident beam normal to the opposite face of the substrate. 

As we proceed, we must take care to adequately quantify the real entities of energy and 
probability as incidence— )>output and output— T-resultant transitions occur. Ordinary incident 
beams generated by conventional sources may be assigned 

P^ = E, (17) 

in arbitrary units which gives an incident beam occupation value 

a = 1 (18) 

as a calculation convenience. The quantification at each transition requires either a verifi- 
cation that Q is maintained or an individual assessment of energy and of probability in that 
transition. 

For any individual slit, the output energy quanta and the output probability necessarily 
remain in constant proportionality relative to that of the incident beam since that single 
slit equally samples both quantities from the incident beam. Indeed, a deviation from this 
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FIG. 2: Emergent output diffractive envelopes from N irradiated individual slits fully determine 
collective output probability (and energy) prior to the formation of resultant orders by interference 
as those envelopes intersect. 

proportionality would constitute a violation of PIQM duality for a mechanism consisting 
of a single slit. Accordingly, from the perspective of LRQM as well as PIQM, an incident 
i7j = 1 is maintained on the output of any single slit as well as collectively over all slit 
outputs, i.e. r^o = 1. See Fig. 2. 

Moreover, with an incident beam of some particular Gaussian diameter D, the (trans- 
mitted) output probabiUty and output energy quanta are both invariant as w and p are 
proportionately varied since 

w/p = a (19) 

presents a constant fractional cross section for transferring probability and energy quanta 
from incidence to output, i.e. the grating's transmission factor. Essentially, proportionate 
increases in w and p results in increases for individual slit outputs of both probability 
and energy quanta but, concurrently, the number of incident-irradiated slits yielding those 
outputs is proportionately decreased. Then, beginning with an ordinary incident beam, 
equal sampling gives the output quantities 

Po;a = E,,, = 1, (20) 

where both are specific to some particular value of a = w/p but independent of w and, as 
a matter of convenience, in arbitrary units we assign unit value to both in Eq. (20). The 
output occupation value VLq is constructed from the ratio of Po-^a and Eo-^a but, because of 
equal sampling for those quantities for any o", their ratios 

Vlo = VLi = l (21) 
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are independent of a. 

The transition to the resultant diffraction orders requires that we examine quantities 
such as wave amphtude and intensity for the emergent output diffraction envelopes. In this 
regard it is most expedient to begin with the classic phasor construction for the output 
wave intensity from a single slit Wos from Kirchhoff diffraction theory in the Fraunhofer 
approximation 

• 2 

Wos{oi) = Wosiof-^ = Wos{0)smc'' a (22) 



where 



a = ^ sin ^ (23) 

A 



for a wavelength A. Eq. (22) above must be used with care as we proceed since the intensity 
Wos{oi) is expressed in "a-space" rather than over the physical 9 angular space azimuthal 
to the slit. Moreover, this phasor construction relates entirely to the slit output intensity 
distribution, giving Wos{(y) relative to a centroid value Wos{0). Effectively then, the Eq. (22) 
intensity is unsealed with respect to the incident beam intensity attenuation in transiting 
the slit of width w. We will return to this consideration below and show that it is readily 
resolved. In the meantime, as a matter of convenience, we choose arbitrary units with 
WosiO) = 1. 

Our present investigation also does not consider slits in the sub- wavelength range w < X. 
Nevertheless, we will have need to consider slit widths in the neighborhood oi w > X where 
classic references such as that of Elmore and Heald caution us that "the Kirchhoff diffraction 
theory is less accurate, and it is expected that the single-slit diffraction factor will no longer 
give a good description of the envelope." [ist] This caution is quite correct. However, we will 
return to this point and show that such deviations are not a factor in the essential criteria 
for duality violation. Accordingly, we proceed here with the use of the sinc^a envelope in 
the interests of illustrating duality violation using an explicit analytical function for that 
envelope and, later in this section, we show that the essential criteria for duality violation 
are well satisfied in the present case and identify the general requirements of an arbitrary 
envelope in meeting these criteria. 

From Eq. (22), we can write the expression for the output single slit probability 

/at 
sinc^ a da (24) 

-at 
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FIG. 3: Intensity plots in arbitrary units of the single-slit output diffractive envelope sinc^a scaled 
by a factor N, the number of irradiated slits, to give a smoothly varying intensity envelope in 
a-space representing the collective output probability of the grating, and the diffractive order 
peaks representing the resultant probability. An untypically small A^ = 4 is chosen here for clarity 
of depiction. The integrals of the two intensity functions are essentially equal, demonstrating 
probability conservation for the a = 1/8 approximation of dense sampling. 



where envelope truncation 



at 



TTW 



(25) 



found from Eq. (23) with 6 = 90° provides the appropriate limits of the Eq. (24) integral. 
In analogy to the collective output probability of two beams incident on a beam splitter, 
the collective output probability of the grating is 



Po-Jat) = N 



sine a da 



(26) 



where here, relative to Pos{at), the coefficient A^ is effectively the number of slits irradiated 
by the incident beam. A^ is well determined for a given beam incident beam diameter D, 
a, and at- The Nsinc^a intensity envelope is shown in Fig. 3. The Eq. (24) and Eq. (26) 
probabilities are identified as functions of the integration limit at- From Eq. (25) we see 
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that those probabihties are equivalently functions of w. Fig. 3 shows an example of an at 
for a narrow sht w > X. The truncation of the sinc^a envelope by ±at at this w results 
in a central lobe and partial side lobes as output from each slit with substantial envelope 
distribution from +90° to —90°. Conversely, significantly increasing at results in a wide 
w ^ X, adds multiple, greatly diminished side lobes to each slit output and concurrently 
confines the significant envelope contribution to the forward direction clustered about 0°. 

However, the Eq. (24) probability Po-a{c(t) computed in a-space, is clearly an increasing 
function oiat, or equivalently, oiw. This result is in confiict with the Eq. (20) w-independent 
probability Po-a deduced from physical considerations for gratings of varying w but having 
some specific a. The confiict can be understood by noting that the Eq. (24) dependence on w 
is an artifactual consequence of integrating the unsealed Eq. (22) intensity in a-space rather 
than in physical space over the 6 angle azimuthal to the slit. This can be appreciated by 
examining the limit of large w for which the w-dependence of Po-a{(^t) becomes vanishingly 
small. In this limit, the significant contributions to the integral are entirely confined to small 
a where a = (nw/X) sin6 — )■ -kwO/X and a in this range is linear in 9. Nevertheless, the Eq. 
(26) collective output probability will prove to be extremely useful precisely because of its 
w-dependence in a-space. 

We next turn to the interference-generated resultant probability emerging from the in- 
tersecting individual slit output sine a amplitude envelopes (shown in Fig. 2) that generate 
the output sinc^a intensity envelopes. With the incident beam of some Gaussian diameter 
D spanning a large number N of grating slits, the non-null resultant probability computed 
in a-space is effectively confined to the sum of the individual integrals of the narrow, highly 
directional principal order intensity peaks, and the intermediate subsidiary peaks are van- 
ishingly small. 

The classic Kirchhoff expression for the resultant diffraction order intensity is given by 



2 /sin(A^a/o') 



WM) = Wos{^)stnc'a .\ / (27) 

\ sm(Q;/a) / 

where VFos(O) = 1 is, as before, the single slit intensity at a = set to unity. The quantity 
(sin(A^Q;/o")/sin(a/o"))^ is the "grating factor". In this form, the diffraction peak maxima 
extend above the single slit sinc^a intensity envelope but are proportional to it. 

In the literature, Eq. (27) is generally modified by replacing VFos(O) with the correspond- 
ing collective N slit intensity VFo(O) at a = and inserting a compensating N'"^ factor 
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based upon the justification that the A^ sht intensity at a = 0, IVo(O) = N'^WosiO). That 
modification would yield the familiar depiction of a diffraction envelope with the maxima 
of the diffraction order peaks in contact with that envelope. However, we deliberately omit 
that modification here in the interests of correctly tracking the relative probability in the 
output— )>resultant transition. Output quantities appropriately refer to pre-interference tabu- 
lations of those quantities. If the single slit output intensity is Wos{0) at a = 0, the collective 
output intensity from A^ slits is NWos{0). The same principle is involved with the Eq. (24) 
single slit output probability and the Eq. (26) collective output probability. 

Consequently, the integral of the Eq. (22) single slit output scaled by A^, i.e. NWos{o:), 
is correctly equal to the integrated Eq. (27) as depicted in Fig. 3. The computed equality of 
these integrals for any selected limits ±at is a confirmation that probability is conserved in 
the output— )-resultant transition. Exact integral equality is achieved as o" — )> for which the 
sample density of the output envelope by the resultant orders becomes infinite. However, 
even the small but finite a = w/p = 1/8 selected in Fig. 3 for illustrative clarity, gives 
an approximation of "dense sampling". In this approximation, integral equality is closely 
achieved as a function of the integration limits ±«f with only minor perturbations as those 
limits pass between discrete resultant orders such as the depicted at displaced from the 
n = 12 order. 

Again, for purposes of illustrative clarity in Fig. 3, an A^ = 4 has been selected, which 
is lower than that for typical experimental conditions by about two orders of magnitude, so 
that the envelope and the diffraction order peaks can both be represented on the same scale. 
Low intensity secondary diffraction order peaks appearing near the bases of the principal 
diffraction order peaks in the figure diminish to negligible values for experimentally realistic 
large A^. 

The integral of the diffraction envelope, which gives the collective output probability 
in a-space, has limits ±at where the envelope truncates at the physical azimuthal angles 
6 = ±90° relative to the grating normal. These limits are a function of w hj at = wn/X. 
The narrow resultant peaks (principal maxima) arise from interference of the intersecting 
single-slit diffraction envelope amplitudes. The peaks occur at a = aj = jna which is jtt/8 
in the figure. The resultant probability can be expressed as a summation in a-space with 
limits ±a;„ where |q;„| < \at\. For small w, the grating is "fine", the limits in a-space are 
low and w > X. For j = 12, 0:12 = 37r/2 and a grating with a truncation value marginally 
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FIG. 4: A detail of a single subinterval of Fig. 3 with A^ = 4 and a = 1/8 showing some arbitrary 



:th 



peak, the associated single slit output intensity envelope and the collective (four) slit output 



intensity envelope with its associated Riemann approximation. The width of the subinterval is 
Act = TTO" and the width of the peak is Actp^ = 7ro"/2 for the present A^ = 4. Most critically, the 
figure depicts the probability equivalence of the peak and the collective output intensity envelope 
over the subinterval. 

inclusive of the ±12*^ orders, i.e. at > 0:12, has w > 3A/2. Conversely, for large w (not 
shown in the figure), the grating is "coarse", the limits in a-space are high and w ^ \. 

Probability conservation in the output— )-resultant transition can also be assessed for indi- 
vidual resultant peaks. From Eq. (27), the grating factor separates the interference maxima 
into subintervals 

Aa = TTcr (28) 

and, within those subintervals, the null to null peak base width is 



Aftpfe = 2TX(y/N. 



(29) 



A detail of a subinterval is shown in Fig. 4. The total output probability for this subinterval 
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IS 

Poj-aictj) = N sin(?a da 

Jaj-Aa/2 

= NiTcr sinc^aj (30) 

where the area under the single sht sinc^a curve at aj is approximated by the product of 

Aa = TTO" and the functional value sinc^aj. (The value of a at the j*^ order is aj = jna 

where the applicable a is implicit from the context.) The total subinterval output probability 

is given by that area, Tcasinc^aj, multiplied by the number A^ of contributing slits. We are 

reminded again that this factor giving total output probability is linear in A^ since it relates 

to an accounting of the collective single slit probability outputs prior to the subsequent 

interference that yields the resultant probabilities. Total output probability in the earlier 

example of two beams incident on a beam splitter provides a useful analog. As a result of 

the linear factor A^ in Eq. (30), the collective output intensity envelope is appropriately 

depicted as the single slit output intensity envelope sinc^aj necessarily scaled up by a factor 

A^ as shown in Fig. 4. This scaling still leaves the resultant diffraction peak maxima above 

the collective output envelope but now graphically illustrates the essential conservation of 

probability in the output— j-resultant transition for this single resultant peak. 

Conservation of probability can also be demonstrated quantitatively for this arbitrary 

single j*^ resultant peak shown in Fig. 4. The probability for this peak is found from Eq. 

(27) noting that 

sin (Naj/a) ^^ 

sm{aj/o-) 
at the principal maxima a = aj. Then 



z, , ^ p+^"-'=/' . 2 f Sin {Na/a)\\ 
Prj;cr{ctj) = / stnc a { — r^ — ^—:— ) da 

Jaj-Aapk/2 

l\ (2tXG \ /,,2 • 2 \ 

\N sine aj) 



sin(Q;/o-) 



2/ V A^ 
= Nil a sinc^aj (31) 

where the integral of the j^^ peak is readily evaluated by noting that, because of peak 
symmetry about aj ± 'Ka/2N , the peak occupies one-half of the rectangular area defined by 
the peak null to null base width Aapk = 2TTa/N and magnitude N'^sinc^aj. Those factors 
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resolve to a resultant subinterval probability and an output subinterval probability, both of 
which can be equated to a Riemann-approximated subinterval area, i.e. 

Prj-a{oij) = Poj;aioij) = 71 a (N siuc^ a j) . (32) 

Accordingly, over all a-space we can define the total resultant probability as a Riemann sum 

n 
Pr;aioin) = ^^ Tl a {N SITK? a j) (33) 

j=-n 

and the total output probability as the corresponding definite integral 

/at 
Nsinc^aj da (34) 

-at 

where Aa = ttc — )> da. We have limit equivalence «„ = at in Eqs. (33-34) when at = mra. 
We can temporarily defer consideration of deviations of the continuous-valued at variable 
from the discrete-valued «„ variable since we are currently imposing dense sampling in a- 
space. Dense sampling divides the integration partition into a large number of subintervals 
{n ^ 1). Consequently, at > an results in fractional integration of Eq. (34) into the 
±{n + 1) subintervals not included in the Eq. (33) summation. The contribution of those 
two fractionally integrated subintervals relative to that of the total 2n + 1 subintervals 
partition is vanishingly small. 

As a formal matter, in the dense sampling limit o" — > 0, the partition subinterval Aa — )■ 
which gives n — )■ oo and 

" pat 

lim 2. ^'^ sinc^ aj = / sinc^ a da (35) 

^ °° „•_ „ J —at 



where the integration limits ia^ are applicable since «„ — )■ a^. The limit condition in Eq. 
(35) is recognized as the Riemann sum equivalency to the definite integral on the right. 
Ultimately, the Riemann sum serves as a mathematical intermediary that demonstrates 
output— )-resultant probability equivalency when a;„ — ?■ at- 

Nevertheless, despite their equivalency, the a-space quantities Po{(^t) and Pr;cr(ttn) are not 
satisfactory expressions of probabilities since an examination of output probability for some 
fixed a predicts w-independence. Both Po{at) and Pj..^„{an) are clearly increasing functions 
of w through their respective w-dependent respective limits at and a„. However, because of 
the relative equality of Po{oit) and Pr-,„{an), the w-dependence of Pr-,„{an) can be removed 
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by normalization with Po{at) which has the same w-dependence. The normahzed resultant 
probability 

En ■ 2 



Pr-A'^t) 



f"' sinc^a^ da 

J -at J 

(36) 



gives a w-independent constant to within a vanishingly small discrepancy arising from at > 
an as discussed above. 

The summation limits ±n are set by the condition |q;„| < \at\. Note that the identification 
of the Eq. (36) normalized expression Pr-cr{(^t) differs from that of the Eq. (33) unnormalized 
expression Pr;cr{<^n) only with regard to the a-space limits. The summation in Pr-a{<yt) has 
the same ±n limits as that of Pr-^a{<^n) but the more distal integration limits \at\ > \an\ of 
the Eq. (34) normalization integral sets at as the defining functional variable of Eq. (36). 
Because of dense sampling, Eq. (36) is very nearly a constant and, as such, is effectively 
independent of at- However, as a formality we retain at as an apparent functional dependent. 
(Since at = irw/X, cti-invariance is equivalent to w-invariance for Eq. (36) where A is a 
constant.) 

From a purely pragmatic viewpoint, the normalization in Eq. (36) expeditiously achieves 
the required ty- invar iance of Pr-a{<yt) for the class of gratings with dense sampling. More 
significantly, however, the Eq. (36) normalization of relative resultant probability by rel- 
ative output probability (distinct from PIQM normalization. Section I) is fundamental in 
quantifying any probability non-conservation in the output— T-resultant transition. For sepa- 
rate real entities of energy quanta and probability, the constant proportionality required by 
PIQM duality is no longer a constraint. In LRQM, deviations from duality may occur as 
non-conservation of probability manifested as a change in the ratio of the resultant probabil- 
ity relative to the output probability. Therefore, for a mechanism that potentially achieves 
duality violation by probability non-conservation in some transition process such as grating 
output— )-resultant, it is critical to represent that resultant probability relative to the base 
output probability in order to express that non-conservation. 

Non-conservation of probability is, however, not in evidence as o" ^ provides dense 
sampling of the outputs by the resultants. Then the Eq. (36) normalized form Pr-cr{<yt) 
rectifies the artifactual increase of the Eq. (33) Pr-a{(^n) with w and provides the physical 
probability as a function of at, albeit a trivially constant unit value. 
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In analogy to the beam splitter with two incident beams, we again have conservation 
of probability in the output— )■ resultant transition and we also have spatial regions of null 
resultant probability and regions of "amplified" resultant probability. Similarly, the mani- 
festation of these resultant probabilities does not alter the reality of the output probability 
wave structures. The grating, however, differs from the beam splitter in that the origin of 
the resultants is physically displaced from the grating surface but still in the near field where 
the emergent expanding individual slit output envelopes begin to intersect and interfere as 
shown in Fig. 2. 

To complete the analysis of the grating with dense sampling, we need to also assess energy 
in the output— 7>resultant transition. As the expanding output probability envelopes interfere, 
the resident energy quanta transfer without loss onto the forming resultant diffraction order 
beams. The process is purely an energy quanta transfer. No wave entity is transferred in 
a process analogous to that of energy quanta transfer at a beam splitter. The consequent 
continued proportionality of total energy and probability in the output— )-resultant transition 
yields a constant Q throughout, 

a = f = fio = §^ = a = |^ = l. (37) 

li 1q If 

We note that in the Eq. (36) Pr;o-(ttt) expression, the property of dense sampling is critical 
to maintaining an 1]^ = 1 independent of w. As at (= rrw/X) is increased, the normaliza- 
tion integral monotonically increases whereas the summation abruptly increases with two 
additional final terms ±n' = ± |n + 1| as a^ reaches a value corresponding to a new j = n' 
order. Nevertheless, for dense sampling, the additional ±n' terms are incrementally small 
relative to the existing sum of the 2n + 1 terms.. 

B. Gratings with sparse sampling 

The dependence upon dense sampling for the maintenance of i7 ?^ 1 suggests that the 
converse, i.e. "sparse" sampling, constitutes a potential mechanism for achieving significant 
deviations of Q from unity which is equivalently a violation of PIQM duality. Specifically, 
we examine a = 0.5, characteristic of Ronchi rulings. This ratio a is notable in that the 
side lobes of the output probability are each entirely expressed by the single, odd order that 
bifurcates those lobes in a-space as shown in Fig. 5. This is a significant fortuitous property 
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FIG. 5: A o" = 0.5 sparse sampling analog to the dense sampling approximation shown in Fig. 3. 
Odd j peaks for j > 3 have probabilities that fully express the probability of the collective output 
lobe that they respectively bifurcate. 

of a (7 = 0.5 grating attributable to the non-expression by the null-valued even orders on 
either end of each of those lobes. It reflects a corollary of Eq. (36) that it is necessarily the 
non-null resultants that express the outputs. The probability for a = 0.5, 



-Pr;0.5("t) 



[vr/2) e;= 



smc a 



J 



(38) 



f"* sinc'^ a da 

J -at 

is merely a special case of the Eq. (36) probability where we again consider w > X (flne 
gratings) to w ^ X (coarse gratings). Similarly, 

(11/2) sinc^ aj 



Prjfl.b{0it) '■ rat ■ 2 J 

\ smc^ a da 

J —at 

for the j^^ order probability. 

However, for the sparse sampling associated with a = 0.5, 



(39) 



a,- 



Jvr/2 



(40) 



and the Eq. (38) probability, unlike the Eq. (36) probability with cr — )■ 0, is no longer a 
constant of the functional variable at- The deviations from constancy occur as at approaches 
the neighborhood of the odd orders. Eq. (38) is plotted in Fig. 6 in a^-space from a^ = tt to 
371. For at = vr, truncation occurs at the ±2"^^ orders, representing a relatively flne grating, 
whereas at = 3n corresponds to truncation at the ih6*^ orders and the grating is modestly 
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FIG. 6: Resultant probability Prfi.bi^^t) plotted here for positive at since functional evaluation 
of the probability at limits at = ±wit/X is understood. The excursions from unity represent 
measurable non-conservation of probability. These excursions are notable for a grating that has a 
w for which at is in the neighborhood of 0^=3. 



coarser. The horizontal axis on the figure shows only positive values since symmetry ensures 
that any choice of at automatically sets integration limits of ±at and concurrently sets 
summation limits ±n in the evaluation of Pr;o.5{c(t) ■ 

An examination of Fig. 6 notably shows a probability edge transition of ~ 5% for at in the 
neighborhood of the ±3'''^ orders. The transitions at the odd orders diminish in magnitude 
as at increases further, ultimately yielding fiat-line constancy of Prfi.nio^t) for very large at 
(coarse) gratings. 

The functional structure of the Fig. 6 probability is, of course, simply a consequence 
of the sparse sampling of the a = 0.5 gratings. As at at aj=2 (±2"'^ order) is increased, 
the Prfi.bio^t) summation is constant whereas the normalization integral increases resulting 
in a maximum ~ 2.5% decrease in the limit as at approaches Q;j=3 from the left. At aj=3, 
-Pr;0.5(<^3) discoutinuously increases by ~ 5% with the inclusion of the j = ±3^'^ terms in the 
summation. As at is increased toward aj=4, the summation with end terms at j = ±n = ±3, 
is again constant while the normalization integral continues to increase resulting in a return 
to the probability value at j = 2. This functional behavior repeats with increasing at 
but with progressively diminishing discontinuities at the odd j. (Since the resultant beams 
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FIG. 7: Theoretical occupation value 0,r-o.5{at) = 0,Qijfyj\i ^^^ for a = 0.5. With energy E set 
to unit value, this function is the inverse of that in Fig. 6. The regions above the unity level 
identify grating w-values that give enrichment whereas those regions below unity correspondingly 
are associated with depletion. The three experimentally measured Q's are plotted for comparison 
on this theoretical curve at their respective j-equivalent values of truncation: ^^^(2.63) ex ^^ ^) 
f^G(3.i6) ex as T, and $1g{3.94) ex as •. 

are physically narrow but finite in Gaussian diameter, Pr;o.5{oit) only approximates a true 
mathematical discontinuity at odd j.) The convergence oi Pr-o,5{at) to a constant, unit value 
in the limit as a^ — )> oo verifies that the Eq. (38) probability is properly scaled. 

Then upon inspection of Fig. 6, we conclude that probability is annihilated in the 
output— )-resultant transition for a = 0.5 gratings with at < aj at odd j. Conversely, 
probability is created in the output— T-resultant transition for a = 0.5 gratings with at > aj 
at odd j. 

As in the above case of at and aj, we have applied a convention in which truncation- 
related quantities are associated with a resultant j^^ diffraction order. It will frequently be 
convenient to return to this convention. For our purposes here, this association is made 
with a positive-valued j but by symmetry the association also applies to the corresponding 
negative- valued j (with, of course, a reversal in the above inequalities). 

We proceed with an examination of the Eq. (38) probability by computing the associated 
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occupation value for this class of gratings, 

-Er;0.5("t) 



^r;0.5(ckty 

f * sinc^ a da 

(41) 



(7r/2) EJ=_„ sinc'^ aj 

since the output energy, which can be set to unity, is conserved in the transition to resultant 
energy, Eo = I = -Er;0.5- Functionally, Qr;0.5{<^t) plotted in Fig. 7, is the inverse of the Fig. 
6 probability -Pr;0.5(«t)- ^r;0.5(«t) then also converges to unity as at — )• oo. For a = 0.5 
gratings with at < aj at odd j, we have enrichment of the resultant diffraction orders. 
Conversely, for a = 0.5 gratings with at > aj at odd j, those orders are depleted. 

It is instructive at this juncture to detail the respective diffraction order parameters P, 
E, and Q for a pair of Ronchi gratings identified as G{3±) that have truncation on either 
side of the significant ±3'''^ orders. Care should be taken to avoid confusion of these two ± 
designations. For ^(3— ), the positive- valued a-space truncation is denoted as «3_ = «* ~ 
as, which marginally excludes the formation of the ±3'''^ order resultant probability channels 
by locating +at at the "interior" (left) null at the base of the +3'^'^ order peak. By symmetry, 
the negative-valued —at is situated at the interior (right) null of the —3'^'^ order peak. 
Similarly, for ^(3+) the positive-valued a-space truncation is denoted as 0:3+ = a* ~ as, 
which marginally allows the formation of the ±3*"^ order resultant probability channels at 
threshold by situating ±at at the respective exterior nulls of the ±3'"'^ order peaks. For a 
given operating wavelength A, the two Ronchi gratings G'(3±) are characterized entirely by 
their respective slit widths Ws±- From Eqs. (25) and (40) at the +3'"'^ order threshold the 
requisite width is ^3 = 3A/2. Therefore the gratings G(3— ) and G(3+) have respective 
widths W3_ ^ 3A/2 and W3+ ^ 3A/2. 

Note that in the above quantities, we have extended the convention of referencing trunca- 
tion to a j'^^ order by use of j— to indicate marginal truncation exclusion of the ±j*^ orders 
and j+ to indicate marginal truncation inclusion of the ±j*^ orders. 

The P, E, and Q parameters for these two gratings G'(3±) are shown in Fig. 8 where 
the collective and the individual j*^ probabilities are computed from Eqs. (38) and (39), 
respectively. The output energy Eo = 1 is conserved giving a total of Er = 1 distributed 
onto the propagating orders as Ej.j. Consistent with probability as a relative quantity, these 
Ej.j are found from the ratio of the Eq. (39) individual j*^ order resultant probability 

25 



E 


=0 
^0.552 




J-0 
E-0.526 




P=0.539 




P==0.539 




\/ 


E-0223 
P=0.21S 


\/ 


E=0213 
P=0,2I8 




/ 


J=-3^ 




4_ 


^ 


Gt3-) 
enriched 


depleted 




j=+3 

E=0-023 

P={).024 



FIG. 8: Tabulation of energy and probability values on the resultant propagating orders for gratings 
G(3— ) and G(3+) where the total energy is unit valued. By symmetry, values on ±j are the same 
for both gratings. 

Prj-o.5{<yt) and the Eq. (38) collective resultant probability Pr-o.5{<yt)- The Fig. 8 example 
illustrates the essential features of probability non-conservation and energy distribution in 
the neighborhood of the ±3^^^ order threshold probability discontinuity. 

We note that since this discontinuity alters the proportionality of the wave-like probability 
and the particle-like energy, the phenomenon can be expressed as a "duality modulation". 
In the Fig. 8 example, the G{3—) grating with ±q;3_ truncation points provides a +2.5% 
duality modulation upon comparing the ^(3— ) occupation value Qg{3-) to the ordinary 
value Q = 1. Similarly, G{3+) provides a —2.5% duality modulation. 

Most generally, this phenomenon of duality modulation is potentially applicable to any 
(fine) grating that yields a small number of resultant orders for which one or two can 
be located near or just beyond the plane of the grating "at threshold". In other words, 
the precise shape of the output diffraction envelope is not a critical factor. This can be 
understood by the relationship of the diffraction envelope and the interference factor that 
modulates that envelope. Within any single subinterval, the j^^ peak arises from the same 
interference factor, identified earlier in Eq. (31), modulating a value f{aj) of an output 
diffraction envelope f{a) in place of sinc^a. Similarly, in the present context of the regular 
Ronchi grating where w > X, the departure from the Fraunhofer approximation w ^ X gives 
a single sht output intensity envelope that is not precisely expressed by sinc^a. However, 
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orders near threshold for gratings with sparse samphng are still predicted to give probability 
non- conservat ion . 

It may be appreciated at this point that the existence of expressed resultant probability 
discontinuities for fine, sparsely sampled gratings could have been deduced directly and 
quite succinctly from basic LRQM principles. The essential rationale for such a deduction 
is derived from the observation that, for a set of gratings with a common sparse sampling, 
e.g. a ~ 0.5, all such gratings have the same output probability, but we can always identify 
specific pairs of these gratings with nearly identical slit widths u'j+ and Wj^ for which Wj+ 
marginally admits a significant resultant order near threshold and Wj„ marginally blocks that 
order. For these two gratings considered in succession, the collective resultant probability 
must exhibit an abrupt decrease in the near infinitesimal Wj+ — )■ Wj^ transition. That 
abrupt decrease does not occur for the smoothly varying output probability. Similarly, the 
output energy is also unaltered by that transition. Therefore, in the absence of an energy- 
dissipative mechanism, the distribution of the output energy onto the resultants yields a 
relative enrichment of the Wj+ resultants with respect to the Wj- resultants. However, 
without the basis developed in the previous subsection for densely sampled gratings, the 
assessment of occupancy in threshold transitions is limited to relative changes in occupancy. 
Moreover, the occupancy is not characterized outside the neighborhood of the Wj threshold 
value. 

Accordingly, we have proceeded here in a less succinct manner by first examining elemen- 
tary beam splitter configurations to elucidate the basic principles of probability conservation 
from the perspective of LRQM. We then examined resultant probability associated with grat- 
ings having dense sampling in order to develop a generalized functional expression applicable 
to all gratings that characterizes resultant probability over the full a-space and not just in 
the neighborhood of threshold transitions before proceeding to an examination of gratings 
having sparse sampling. 

It is of some peripheral interest that the spatial physical discreteness of resultant diffrac- 
tion orders of a grating is the essential property that produces duality modulation in the 
present example. We recall that it was also spatial physical discreteness that yielded dual- 
ity modulation for the case of a discrete photon beam incident on a beam splitter. For a 
grating, it is the discrete increment of probability (on a highly directional diffraction order) 
that modulates duality whereas for the beam splitter, it is the discrete increment of energy 
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that produces this modulation. 

We conclude this section on the theoretical basis for duality violation with a brief but rel- 
evant examination of grating anomalies. This examination is necessitated by some apparent 
similarities of those anomalies to the phenomena considered here. 

Grating anomalies, since their discovery by Wood in 1902, have continued to be a subject 
of intensive experimental and theoretical investigation, [ij] These anomalies are deviations, 
often abrupt, in the diffraction order irradiances as a function of a parameter such as wave- 
length. These deviations were initially designated as anomalies since they departed from 
predictions of earlier classical principles thought to fully characterize grating phenomena. 
Their continued characterization as "anomalies" is now generally regarded as a misnomer 
since sophisticated theoretical electromagnetic analyses of these phenomena have substan- 
tially approximated the experimental observations. 

Our particular interest here concerns "Rayleigh" anomalies that have a theoretical ba- 
sis postulated by Lord Rayleigh several years after their discovery by Wood.|l5| Rayleigh 
anomalies are characterized by abrupt increases in the irradiance of propagating orders as 
one of those orders "at threshold" near the grating plane is extinguished by an incremental 
wavelength increase. The essential theoretical basis postulated by Lord Rayleigh consists 
of coherent (photon) scattering of the threshold order on a grating's periodic structures. 
By interference, the scattered threshold order (inclusive of its irradiance) is coherently re- 
distributed onto the remaining propagating orders in proportion to the relative intensities of 
those orders. This coherent scattering process is fully consistent with PIQM which was devel- 
oped two decades later. The postulated theoretical basis transfers photon energy quanta as 
well as photon wave probabilities to the remaining propagating orders thereby maintaining 
duality. 

We seek here an alternative demonstration in LRQM of Rayleigh-like anomalies in the 
absence of the coherent scattering process postulated by Lord Rayleigh. With a = 0.5, the 
slit width w is selected as the dependent parameter consistent with our previous analysis. 
From Eq. (25), for which w and A are inversely related, LRQM predicts abrupt irradiance 
increases as orders reach threshold with progressive incremental decreases in w. 

In the open interval between any successive j = n and j + 1, each of the individual 2j + 1 
propagating orders exhibits a proportionately decreasing probability across that interval as 
shown in Eq. (39) and Fig. 6. (As in the Sec. 3 analysis, we consider the positive orders, 
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FIG. 9: order energy as a function of j(ty)-equivalent truncation showing characteristic edge 
drops in propagating energy (or irradiance) associated with Rayleigh anomahes at threshold values 
of significant orders. Compare to Fig. 8. 

but the corresponding negative orders are inclusive by symmetry for normal incidence.) The 
resultant probability of any single order relative to the total resultant probability Eq. (38) 
is constant for any at within any such open interval. This fractional resultant probability 
for the O*'' order 
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(42) 



is found from the ratio of Eq. (39) with j = and Eq. (38). The independence of this 
ratio with respect to the a-space location at is emphasized by that ratio's equivalence to 
evaluation at the initial a-space location a„ of that interval, i.e. i?o(ctra)- Essentially, the 
Qth Qffjgr'g fractional share of the total resultant probability is calculated with respect to 
the probability of all propagating orders as a function of increasing at- This fractional share 
is a step function, constant for some total number of propagating orders 2n + 1 and then 
discontinuously dropping as at reaches the next order at some j = n + 1. The 0^^ order's 
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fractional share, calculated against the total probability of the 2(n + 1) + 1 propagating 
orders, drops because of the additional probability on the ± |n + 1| orders. Fig. 9 shows a 
plot of the 0*^ order fractional share /?o(<^t)- 

In LRQM, energy quanta entering a probability field equilibrate onto that field in pro- 
portion to relative probabilities. Since the total output energy quanta in that field are 
independent of at-, the O*'^ order's fractional share of total resultant probability i?o(ttt) also 
gives the fraction of the total output energy Eo quanta equilibrating to the O*'' order. Then 
the energy on the O*'' order is 

E,o = E,R^{at) (43) 

and with ii^o = 1 in arbitrary units, -E^o is simply represented by the Fig. 9 RQ{at) step 
function. This function predicts the abrupt step-wise increases associated with Rayleigh 
anomalies both with regard to their locations at threshold at and their magnitudes relative 
to the irradiance of the paired extinguished threshold orders. 

However, unlike the (PIQM-consistent) coherent scattering process postulated by Lord 
Rayleigh, the LRQM-consistent basis described earlier in this section also predicts a du- 
ality violation as quanta transfer from an extinguished threshold order to the remaining 
propagating orders without concurrently transferring wave intensity to those orders. 

We stress that this finding does not preclude the existence of PIQM-consistent postulated 
coherent scattering as the operant mechanism for Rayleigh anomalies in other gratings and 
beam configurations, but it does expand the potential mechanisms that may be the cause of 
observed anomalies to include those that are not consistent with PIQM. In this regard we 
note that Lord Rayleigh based his theoretical analysis on experimental evidence from refiec- 
tive gratings having deep, pronounced periodic groove structures which would be expected 
to exhibit significant scattering of a threshold order into the field of the remaining prop- 
agating orders. Conversely, the transmissive gratings we consider here consist of periodic 
fiat opaque bands of negligible thickness that would not provide a comparable scattering 
mechanism. 
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FIG. 10: Experimental apparatus configuration showing a potentially duality-modulated beam ^q 
following passage through a particular grating G. ^g is coupled to an ordinary beam <l>/j along a 
coupling path extending from a beam splitter BS. Coupling occurs with beam blocker B^i shifted 
to transmit ^r. 

IV. EXPERIMENT 

A. Apparatus and beam parameters 

In this subsection we begin with a description of the experimental configuration shown 
in Fig. 10 before proceeding in subsequent subsections to the methodology for assessing 
duality violation on that configuration. 

A HeNe laser generates a horizontally linearly polarized beam $ of several milliwatts at 
633 nm. Beam $ traverses a variable attenuator Att and an optical beam chopper wheel 
Ch. $ is normally incident on a grating G. The grating is one of several Ronchi rulings with 
various grating slit widths. The intervening opaque bands of the gratings defining those slits 
are 150 nm-thick refiective chromium deposited on a glass substrate. The particular grating 
under study is mounted with the ruling on the exit face and the bands vertically oriented. 
The 0*^ order diffraction beam identified as ^q is incident on a 50:50 beam splitter BS. An 
independent, HeNe laser generates a horizontally linearly polarized beam ^/j initially several 
milliwatts in power. (^^ passes a retractable beam blocker B^i and enters a beam expander, 
L 1 (/ = +100 mm) and L2 (/ = +200 mm), before forming a beamspot concentric with 
that of <I>G on the beam splitter BS as shown in the Fig. 11 detail. This concentricity is a 
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FIG. 11: Detail of Fig. 10 coupling path showing the substantial separation of $r from ^g at the 
path terminus by convergence of the former onto a disk mask. 

critical alignment for the apparatus. As a practical matter, the experimental configuration 
includes a number of beam directors not shown in Fig. 10 that facilitate beam alignment 
and folding of optical paths. 

As we shall see in the subsequent subsections, measurement of duality violation imposes 
some general criteria on the beam and apparatus parameters. However, in the interests 
of facilitating replication of the present experiment, we provide specific parameters here in 
greater detail than that necessitated by the general criteria. 

The (Gaussian) diameter of $g at BS is ~ 2mm by natural divergence from the source 
laser and whereas the corresponding diameter of $ij at BS is ~ Am,m, as a result of LI and 
L2. The beam components exiting BS utilized here are the transmitted component of ^g 
and the reflected component of ^r. The orientation of BS is adjusted to coaxially align the 
^R beam spot to the ^g beam spot at the terminus of a 2000mm "coupling path". This 
critical, second coaxial beamspot alignment effectively coaxially aligns the $/? and ^g beams 
on the coupling path. The variable attenuator Att significantly reduces the beam power of 
$ such that ^G along the coupling path is ~ 6fiW. The corresponding beam power of $k 
along the coupling path is approximately two orders of magnitude higher at ~ 600;uW^. ^g 
expands to a Gaussian diameter of ~ 4.5mm at the terminus of the coupling path by natural 
divergence from the emitting laser whereas $« converges to < 1.7mm by setting the relative 
spacing of lenses LI and L2. 

An assembly of a beam blocker Bm, an iris diaphragm Ir, and a photodiode detector 
Det is located at the coupling path terminus. Bm consists of a 1.7mm opaque disk mask 
mounted on a transparent glass substrate. The iris is set to a 3.3mm diameter. Beam 
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directors on the coupling path (not shown in Figs. 10 and 11) are used to provide concentric 
ahgnment of $g and $/? with Bm and Ir. These settings colhmate the ~ 6fiW $g beam on 
the couphng path to an annular beam of ~ 2A^W incident on the adjacent detector Det. 
Of the remaining power, ~ l.b^W (25%) is blocked by the disk mask B„i and ~ 2.1fiW 
(35%) is blocked by the diaphragm of Ir. The ~ GOOfiW power of $/? on the coupling path 
is blocked almost entirely by Bm alone leaving a residual ~ 2.5fiW incident on the detector. 
The function of the iris setting is to simultaneously restrict the radial sampling of $g to an 
annular region closely coupled to $/j while concurrently providing adequate ^g power for 
detector measurement. 

B. Observation of duality violation 

We briefly digress in this subsection from the experimental apparatus and consider the 
theoretical basis for observing duality violation. As we discussed in Section III, orders 
passing to threshold for particular gratings are associated with abrupt redistributions of 
irradiance onto the remaining propagating orders for Rayleigh anomalies. A similar redistri- 
bution is predicted here for particular Ronchi gratings but, additionally, with the predicted 
phenomenon of duality violation on those remaining propagating orders. Consequently, for 
these Ronchi gratings the duality state of the propagating orders must be measured in order 
to establish whether those orders are ordinary, as would be expected for Rayleigh anomalies, 
or are in fact duality modulated. The experimental design used to achieve this measurement 
requires a transient coupling between a resultant presumptively duality- violating beam from 
the grating (a propagating order) and an independent ordinary beam. [6 1 

This design is analogous to the intersecting of two independent beams as used in numerous 
investigations 16|] to experimentally assess duality violation by determining the presence or 



absence of interference between those beams, i.e. to provide a test of PIQM. A review of 



these investigations is given by Paul. 17(| An observation of interference would seemingly 



violate Dirac's dictum that a photon (in PIQM) can interfere only with itself. |18| 

The outcomes of these numerous investigations are conclusive demonstrations that inter- 
ference does occur in apparent contradiction to PIQM. However, in a theoretical analysis of 
this phenomenon, Mandel makes the critical argument that for any given photon measured 
in the interference we do not know on which beam that photon had initially resided. [l9| 
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Because of that lack of knowledge, each photon is treated in Mandel's analysis as inter- 
fering with itself. Consequently, interference in the intersection of independent beams is 
consistent with PIQM and does not provide a test of that interpretation. Concurrently, that 
interference is trivially consistent with LRQM. 

In a variant of those independent beam investigations, we prepare one of the beams by 
transmission through a particular grating where we have some basis that a prepared beam 
<I>G Kiay specifically violate PIQM duality, i.e. the beam is in a depleted or enriched state 
from the perspective of LRQM but is necessarily ordinary for PIQM. Spatially transient 
equilibration of that prepared beam with a second, ordinary beam ^^ by mutual interference 
over the coupling path should then yield a net transfer of energy quanta for LRQM but not 
for PIQM. That net energy transfer relative to $g is readily measurable in the cw regime 
by detecting disparate beam powers on the sampled $0 with and without ^r present on 
the coupling path. 

C. Gratings 

The particular gratings used in the experiment are based upon the theoretical analysis 
presented in Section III. In the apparatus shown in Fig. 10, the horizontally linearly polarized 
(A = 633nm) $ is normally incident on one of three Ronchi transmission gratings G with 
respective slit widths w that are theoretically predicted to generate resultant orders that are 
respectively enriched, depleted and ordinary. For whichever of the three gratings is installed 
in the Fig. 10 apparatus, the lines on the exit face of that particular grating are vertically 
oriented thereby providing TM (S) polarization with respect to G in the usual classical 
configuration for observing grating anomalies. 

From Eq. (25), for a specified incident wavelength, the slit width uniquely determines the 
truncation point at- From the theoretical analysis in Section III, the most relevant attribute 
of a particular Ronchi grating is its at truncation point with respect to the location of 
some particular j*^ order. The three Ronchi gratings used in this investigation are uniquely 
characterized by the respective slit widths w = 833 nm, 1000 nm, and 1250 nm. From Eqs. 
(25) and (40), j{w) = 2w/\ gives a continuum of w-dependent truncation points relative to 
the j^^ integer diffraction orders that is 

j{w) = 2.63, 3.16, and 3.94 (44) 
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for the three respective sht widths. In the present investigation, A = 633 nm is invariant. 
Accordingly, the individual gratings are also uniquely characterized by the Eq. (44) con- 
tinuum j(tf;)-equivalent truncation points. A Ronchi grating of some arbitrary slit width w 
is denoted as G[j{w)] or simply G. Conversely, a grating identified with a numerical j{w)- 
equivalent truncation point identifies a particular grating with an implicitly expressed slit 
width. This use of j(u') continues the convention in which variables are most instructively 
identified by the critical j-equivalent truncation value. 

In this convention we are reminded that a grating such as G(3) has truncations precisely 
bisecting the ±3'''^ orders whereas a grating G(3— ) has truncations marginally exclusive of 
the ±3'''^ orders i.e. truncations at the respective interior nulls of the ±3^^^ peaks. Extending 
this to the j{w) continuum, a particular grating G(2.63) has an implied slit width w = 833 
nm. The truncation point at j = 2.63 is exclusive of the ±3*"^ orders but only approximately 
satisfies the marginal exclusion of those orders specified by j = 3—. Similarly, a grating 
G'(3.16) implies w = 1000 nm and includes the ±3'''^ orders, but the truncation point at 
j = 3.16 is not marginal as it is for j = 3+. Despite the lack of truncations marginally 
close to j = 3 for G(2.63) and for G'(3.16), these gratings are nevertheless predicted to yield 
orders that are respectively significantly enriched and depleted. 

The final grating, G(3.94) with w = 1250 nm, is included to provide for a control exper- 
iment since the j-equivalent truncation point very closely matches the null orders at j = 4. 
See Fig. 7. Gratings with a j-equivalent truncation in the immediate neighborhood of j = 4 
are predicted to provide orders that are ordinary. Consequently, G(3.94) orders are expected 
to have no significant net energy transfer from coupling and should then exhibit no duality 
violation. 

D. Calculation methods 

Data are acquired with one of the three Ronchi gratings in the Fig. 10 position of G. 
Since the Gaussian beam diameters of the coupled $g and $/? are not equal over the coupling 
path as a result of lenses LI and L2, we necessarily undertake the examination of coupling 
phenomena with the extensive variables E and P rather than their respective intensive fiux 
densities irradiance / and wave intensity W^ = |<I>| . We continue the use of ^g and $k as 
general identifiers of the respective beams, but we are reminded that these wave functions 
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in LRQM are exclusive of the energy quanta residing on those beams. 

The basic premise of beam couphng is that a duahty modulated beam equilibrates with 
an ordinary beam by a net transfer of energy quanta that leaves the wave structures of both 
beams unchanged and ideally converges the occupation values toward a common value. 
This convergence objective constitutes criterion (1). For this idealized coupling, a fully 
equilibrated state is achieved as $g and $/? approach the end of the coupling path giving 
the equality 

^Gc = ^Rc- (45) 

We use the added subscript "c" on variables such as Q to denote values at the end of the 
coupling path where equilibration of $g and $/? has potentially altered those values relative 
to their respective values without $g and ^ji simultaneously present on the coupling path 
for those same variables. 

Under a criterion (2), ^ji should ideally serve as an infinite source for a depleted <I>g or 
an infinite sink for an enriched $g in the equilibration process. In the present case with 

satisfied by a respective ratio of ~ 100 : 1 for these respective probabilities, we have an 
approximation of criterion (2) leaving the final equilibrated ^g and $k both as ordinary 
and extending the Eq. (45) Q equality to a unit-valued ordinary value, 

^Gc = ^Rc = 1. (46) 

If $(7 is depleted or enriched as it emerges from a grating G, a net transfer of energy AE 
will occur between ^g and ^^i that changes the initial grating-emergent energy Eq of <I>g 
to 

Egc = Eg±AE (47) 

as the coupling path terminus is approached. A positive-signed AE corresponds to an energy 
gained by ^g in a transfer from ^r where ^g had initially been depleted. Similarly, if $g< 
had initially been enriched, AE is negatively signed. Alternatively, if ^g emerging from G 
is initially ordinary, no net transfer occurs and AE is zero. 

The coupling equilibration of ^g to an ordinary state (if it is not already in an ordinary 
state) provides at the coupling path terminus in our arbitrary units the important result 

Pg = Egc. (48) 
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FIG. 12: Graphical depiction of irradiances at the couphng path terminus showing, in particular, 
the respective ^g ai^d ^r annular beam irradiances Ic-a and Ir^u incident on the detector Det. 



A sampling of this Eqc is readily acquired by the detector. 

Under idealized criterion (3), ^r is entirely excluded from the detector annular sampling 
region. In the present Fig. 11 configuration with the given beam parameters, the ~ 100 : 1 
ratio of ^r probabilities PR-m blocked by the B^n mask and Pr-u incident on the detector 
sampling region approximates criterion (3) as depicted in the Fig. 12 graph of irradiances. 

The Fig. 13 oscilloscope pattern for the detector sampling with the chopper wheel in 
motion is a square wave of the pulsed $g energy with a baseline biased by some steady state 
energy on ^r residually in the annular sampling field as well as any background level. Then 
the oscilloscope peak height measurement of the square wave 



AVbc = k^Egc 



(49) 



gives the beam's post-coupled energy and the ^g probability as well because of Eq. (48) to 
within a multiplicative constant k,. Chopper wheel pulsing of the $g< beam enables a peak 
height measurement of the square wave that automatically separates the detector sampling 
of Egc from any steady-state bias sources. For acquisition with a particular Ronchi grating 
G in position, the chopper wheel transmits and blocks <I>g for equal 5 msec time intervals 
giving 10 msec cycles in generating a square wave pulse train transmitted to the oscilloscope 
from the detector amplifier. 
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FIG. 13: Square-wave pulse outputs from detector amplifier (not to scale) for an initially depleted 
^G- The measured pulse height AV^ with Br blocking ^^ increases to AVqc by equilibration 
transfer of irradiance from <!>/? on the coupling path. 

After AVgc is acquired, $/? is blocked from the coupling path by B^. The detector then 
samples the same annular region of <I>g but now the peak height measurement 



AVn = nE, 



G 



(50) 



provides the beam's grating-emergent energy Eg-, unmodified by coupling, to within the 
same multiplicative constant n. The vital significance of these two detector measurements 
is that their ratio 



AVc 



G 



AVf 



Gc 



Eg 

^G ex 



(51) 



which is the experimentally determined occupation value from a single pair of measurements 
AVgc and AV^. The subscript "ex" has been added to clearly identify this quantity as 
experimentally determined. The final averaged Qg ex reported in the next subsection for 
each of the three gratings studied are computed from multiple Eq. (51) determinations 
which are themselves averages of ~ 100 pulse cycles. 

The theoretical counterpart to those final Qg ex experimental values is the Eq. (41) 
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function 

fir;0.5(ttt) = rir;0.5(7rw/A) 

= ^G[j(w)] th, (52) 

shown in Fig. 7 and re-designated here as Qg[j{w)] th since the wavelength is invariant at 
A = 633 nm, G is uniquely characterized by the continuum value j{w), and the added 
subscript "t/i" explicitly emphasizes the theoretical basis of the function. 

In the foregoing method for calculating Q from experimental measurements, it is necessary 
to cast the equations in terms of the directly detectable beam energies (as derived from 
beam powers) in order to infer the beam probabilities which are themselves not readily 
amenable to direct measurement. However, it is probability creation and annihilation that 
is of fundamental significance to duality violation in LRQM and therefore it is important 
from a theoretical perspective to explicitly express the resultant occupation value entirely 
in terms of probability. 

For the general conditions of an ordinary output 

"o p -L 



and a resultant 

With energy conservation 

in the transition and identifying 

we have the equation 



JTr 



Eo = Er 



Pr = Po± AP, (53) 



n,.= [l± — j (54) 

where a positive sign represents AP probability creation and a negative sign represents AP 
probability annihilation. 



E. Experimental results 

Consistent with the notation in Eq. (52), the experimentally determined Qcijiw)] ex values 
specific to the three gratings are f^G(2.63) ex, ^0(3.16) ex, and f^G(3.94) ex- 
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The three experimental values Qg ex are plotted for comparison on the Fig. 7 theoretically 
predicted Qg[j(w)] th- The experimental values ^G(2m) ex = 1-015 ± 0.003 with a duality 
modulation of +1.5%, ^0(3.16) ex = 0.982 ± 0.003 with a duality modulation of —1.8%, 
and ^0(3.94) ex = 0.998 ± 0.003 with a duality modulation of —0.2% are in good agreement 
with the theoretically predicted VLg[j{w)] th function for LRQM shown in Fig. 7. At the 
most fundamental level, quite independent from the LRQM theoretical basis presented here, 
any experimentally significant VLg^j^^-^-^ ex deviations from unity demonstrate net transfers of 
energy during coupling that are inconsistent with PIQM. 

F. Analysis of biased error 

We conclude this section with an analysis of the three idealized coupling criteria identified 
in subsection D: 

(1) The coupling path has perfect equilibration efficiency. 

(2) $^ serves as an infinite source or sink. 

(3) ^R is totally excluded from the detector sampling region by the mask on Bm- 

As a practical matter these criteria are not fully achieved in the Figs. 10 and 11 apparatus 
and, accompanying the usual statistical dispersion of measured values, there are biased 
(non-random) sources of error present related to these criteria. We show below that these 
sources cause the experimentally measured VLq ex to underestimate the duality modulation, 
i.e. the actual magnitudes of the duality violations are larger than that of the current 
experimentally measured duality modulations. Moreover, we show that the magnitudes of 
the underestimates are not a significant fraction of their respective duality modulations. 
Despite the smallness of these underestimates, we include the biased error analysis here as 
an exercise in completeness and in the interests of identifying how the measurements can 
be optimized. These criteria and calculation of Vt underestimates are examined in greater 
detail in ref. [6J. 

Criterion (1) relates to an idealized coupling of two beams that fully equilibrates them 
to a common Vt. This equilibration is mediated by the physical proximity of the two beams 
and the longitudinal extent of that proximity. Both of these factors are addressed in the 
current experimental configuration by the coaxial propagation of these beams on a 2000 mm 
coupling path and the iris diameter. 
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Matching of the respective beam Gaussian diameters along the couphng path is necessarily 
altered by deliberate convergence of ^r onto the Bm mask in order to provide substantial 
separation of that beam from the detector annular sampling region. Nevertheless, this 
configuration still provides for effective coupling over the 2000 mm path. Longer coupling 
paths have been examined in the interests of potentially improving the coupling efficiency, 
but with no significant change in Qq ex observed, the 2000 mm path was retained. 

Similarly, for the Fig. 11 beam configuration, coupling path efficiency is further optimized 
by maximizing radial equilibration between $g and $/?. This is achieved by confining the 
sampled $g to an annular region most closely coupled to the convergent $ij using a minimal 
iris setting (that still admits sufficient $g beam power for measurement). 

In any case, an incomplete equilibration arising from a non-ideal coupling implies that 
the AE transfer is less than it would be if Qqc — > 1- As a result, the apparent Pq equated to 
the measured Eqc would be underestimated for depletion and overestimated for enrichment 
yielding an underestimated magnitude of duality modulation for both. 

In order to demonstrate the magnitude of the duality modulation underestimate for the 
not fully achieved criteria (2) and (3), it is most useful to begin with idealized experimental 
values ^G(3-) ex,id = 1-025 and ^^(3+) ex,id = 0.975 for which criteria (2) and (3) are achieved 
in principle and calculate the respective apparent values f2G(3±) ex,ap- 

For criterion (2), if $/j has only a finite probability Pr, perfect coupling path equilibration 
from criterion (1) does still provide equalized final occupation values 

^G(3±)c ex = ^Rc ex,id 7^ 1 (55) 

but because of the finite Pr, the initial i^G(3±) ex,id and Vt^ are mutually convergent on a 
non-unit value. For the intensive fi's and extensive P's and E^s, it can readily be shown 
that 

_ Eg + Er 

^l'G{3±)cex - p , p 

= na.x,.dp^^+^np^^ (56) 

where, for the initially ordinary <l>/j, Er = Pr and VIr = l.|6| 

To examine this expression in the context of the present experiment, where Er ^ IOOEq, 
we note that the probabilities are also approximately related by Pr ~ IOOPg despite a small 
inequality of Eq and Pq for a small duality modulation of the initial $g- These energies, 
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which are measurable by detector, may be substituted for the respective probabihties in the 
coefficients Pg/{Pg + Pr) and Pr/{Pg + Pr)- This substitution introduces only a second 
order error in Eq. (56) since the respective P ^ E. Then 
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(58) 

Then, for the finite $/? that does not yield a fully realized criterion (2), Egc = 1.0002Pg is 
treated as Pg and the apparent experimentally measured occupation value 

fiG(3-) ex,ap = ^^ = 1.0002Pg ^ 10002^^^='-) '^"•*"- ^^^^ 

From this result we find that the apparent f2G'(3_) e^.ap is insignificantly smaller than the ide- 
alized ^G{z~) ex,id and the apparent duality modulation ^g{3~) ex,ap — 1 is also insignificantly 
smaller than the idealized value. 

Similarly, for i^G(3+) ex,id the apparent ^2^(3+) ex,ap is insignificantly larger than the ideal- 
ized ^G{3+) ex,id and the magnitude of the apparent duality modulation |fiG(3+) ex,ap — l| is 
also insignificantly smaller than the idealized value. 

Lastly, we consider deviation from criterion (3). In practice, the Gaussian tail of the 
convergent $/? beam spot always has some small but finite fraction of Pr outside the B^ 
mask into the annular sampling region. As a consequence, $/j is not completely separated 
from the annularly sampled ^g- Conversely, that sampling accepts a substantial fraction 
of Pg- The general significance of this incomplete separation relates again to an apparent 
discrepancy in AE transferred to or from $g fo produce an Egc that is presumed to have 
equivalence to Pg- In the equilibration process an equal AE is transferred respectively 
from or to $/?. Upon coupling, some large fraction Eg of AE transfers to or from $g in 
the annular region while concurrently some small but finite fraction Pr of AE transfers 
respectively from or to $^ in the same annular region resulting in a diminished net energy 
transfer in the annular measurement region. 
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This dependency of the occupation value on the portion of ^r in the annular sampling 
region is given by an apparent 

^G(3±) ex,ap = ^ j^, _^ ^ ^ ^ ^ c. a z:^\ (^0) 



FqEg ± (FgAE - FrAE) 



■E 



where ± applies respectively to an initially depleted or enriched $g-|6| To put this error into 
perspective with regard to the apparatus of the present experiment, Fq ~ 0.4 and Fr ^ 0.01. 
For the idealized ^^(3-1-) ea-j^, the magnitude of the duality modulation is 1^2^(3-1-) ex,id ~ l| = 
^E / Eg = 0.025. With these values, an idealized 2.5% magnitude of duality modulation 
for both depleted and enriched $g is again underestimated as an apparent insignificantly 
smaller magnitude. 

We conclude that the close approximations to the criteria (1), (2), and (3) for the beam pa- 
rameters in the present experiment provide apparent VLq ex,ap values that differ insignificantly 
from idealized Qq ex,id values. Consequently, no compensating adjustment of experimental 
results Qq ex in the previous subsection is required with regard to the coupling criteria. 

V. DISCUSSION 

A viable locally real alternative to the probabilistic interpretation of quantum mechanics 
PIQM must necessarily be in agreement with performed experiments and must provide a 
self-consistent theoretical basis for representing quantum mechanical phenomena. 

We hypothesize that a comprehensive locally real representation of an underlying quan- 
tum mechanical formulation, LRQM, can be systematically constructed by providing for 
separable locally real wave-like and particle-like physical entities where the wave-like entity 
is represented on a covariant field of elementary oscillators generally in ground state with 



.0 



relative non-coherence. |^ A photon in this representation consists of a wave packet of these 
ground state oscillators in relative coherence, i.e. the wave-like entity. The particle-like 
entity, i.e. the energy quantum, consists of a raised energy state of one of these oscillators 
in coherent motion. The density of these oscillators in relative coherence at a particular lo- 
cation on the wave packet $ is given by |$| . This density of oscillators in relative coherent 
motion represents a small fraction of the total density of field oscillators at any given point. 
Consistent with Born's rule, I'I'I provides the probability flux density for the location of 
the energy quantum on the wave packet. For ordinary (0 = 1) photons, the integration of 
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|$| over all space (essentially, over the entire wave packet) yields a value that is in strict 
proportion to the unit energy quantum consistent with the PIQM principle of duality. 

We emphasize that these LRQM interpretations, restricted to a simple system such as an 
ordinary discrete photon, are not measurably distinguishable from those of PIQM. Differ- 
entiability arises when the structure of the wave packet is altered in such a way that local 
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probability on a photon is not conserved. For correlated entities, this is shown in ref. [^ to 
be a constraint on the ensemble of quantum states in Hilbert space that results in a proba- 
bility loss on one of the correlated entities. The locally real solution, which is independent 
of Bell's Theorem [5| , is consistent with the underlying quantum mechanical formalism and 
agrees with performed experiments. 

In the present investigation of duality, non-conservation of probability is again the factor 
distinguishing LRQM from PIQM. With respect to duality, the probability non-conservation 
is manifested as an alteration of the local flux density I'I'I of oscillators in relative coher- 
ence. At the simplest level, the wave packet of an ordinary discrete photon traversing a 
beam splitter emerges as a pair of spatially similar wave packets along the transmissive and 
reflective output channels but with relative oscillator coherence densities reduced to T|$| 
and -R|$| , respectively where T is the transmission factor and R is the reflection factor 
of the beam splitter. The energy quantum transfers onto one of these two packets, again 
consistent with Born's rule. Probability is conserved when summed over both output chan- 
nels [T + R = 1) and PIQM retains duality by invoking a non-local probabilistic (non-real) 
interpretation of this phenomenon. From the perspective of LRQM, the two emergent pack- 
ets, only one of which is occupied by the energy quantum, both exhibit physical distinctions 
from the incident photon. However, the testability of these distinctions at the discrete level 
is problematic. This difficulty is circumvented here by going to the continuous wave (cw) 
regime and using a grating system presumptively generating duality-modulated beams since 
even very modest duality modulations are definitively testable in the cw regime. 

The LRQM treatment of probability as a relative quantity can be demonstrated from 
the perspective of the densities of the field oscillators in coherent motion on the resultant 
diffraction beams emergent from a grating near threshold. For example, ~ 2.5% of the 
output probability that forms the resultant probability is annihilated in the near zone of 
grating G{3—). Equivalently, this translates to a ~ 2.5% reduction in the coherent oscillator 
density on each of the ^(3— ) resultant orders. However, as the output energy quanta in 
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the near zone of G{3—) transfer onto the resultants, the proportionate ~ 2.5% reduction 
in the coherent oscillator density on each does not alter the relative distribution of the 
total quanta onto those orders. Indeed, any proportionate reduction or increase in coherent 
oscillator density on a complete set of resultant channels does not change the distribution 
of a given set of quanta transferring onto those channels. Probability, treated as a measure 
of the physical coherent oscillator density, is appreciated as a relative quantity expressing 
the expectation of quanta transferring to a particular channel. |9| 

Conversely, the motivation to compactly fold conservation of the particle-like energy 
quantum in with the distributional probability results in a PIQM profoundly distinctive from 
LRQM. The PIQM normalization of the integrated I'I'I over all space effectively elevates 
the wave-like probability to equivalence with the particle-like energy quantum. This is really 
the fundamental statement of duality. An objective review of the distinctions imposed by 
PIQM duality is given by Rabinowitz. (20[] 

In PIQM energy quanta and probability are then dual manifestations of a single entity, 
the "photon" . However, when superposition states such as those created by a beam splitter 
are considered, PIQM is forced to impose the properties of non-locality and non-reality 
which particularly distinguish that interpretation from LRQM. 

Alternatively, in the comprehensive LRQM that emerges, the field coherence states are 
represented by the wave functions that exhibit probability non-conservation for particular 
quantum phenomena. These wave functions, freed from the constraint of duality-imposed 
renormalization, "complete" the underlying quantum mechanical formalism in the regard 
that the quantum phenomena are then representable as locally real. In this context, Einstein, 
Podolsky, and Rosen had referred to objective properties of a physical system that are 
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represented by parameters which they called "elements of reality". |2| In common usage 
these elements have since, unfortunately, been re-identified as "hidden variables". Ferrero, 
Marshall, and Santos present a compelling argument for the inappropriateness of this re- 
identification. Their argument has application here since the LRQM wave function is shown 
to express objective properties before any measurement is made on it as they prescribe in 



ref. 2l|. Moreover, the objective properties of that wave function can realized by physical 



measurement. The intensity |$| of a beam is disproportionate to that beam's irradiance 
for a non-ordinary beam. That disproportion is measurable from the energy quanta transfer 
that occurs during equilibration coupling of that beam with an ordinary beam. 
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From the perspective of PIQM, for which probabihty is systematically conserved in the 
preparation of the normalized wave function, the variables that would provide locality appear 
to be hidden since phenomena that distinguish PIQM from LRQM are associated with non- 
conservation of probability. 

Ferrero, Marshall, and Santos postulate that "in spite of the spectacular success of quan- 
tum mechanics, it is worthwhile exploring (small) modifications of the formalism in order 



to ensure compatibility with local realism." 2l| This comment was made in the support of 



seeking a locally real theory that naturally obviates entanglement, but necessarily implies 
extension to an encompassing theory that also naturally excludes quantum phenomena that 
invoke non-locality such as duality. The objective posed by Ferrero, Marshall, and Santos 
is intrinsic to LRQM by incorporating the underlying quantum formalism modified only by 
the omission of probability renormalization in transitional processes. The resultant LRQM 
is broadly consistent with PIQM predictions and performed experiments while providing for 
a local reality-based representation. 

VI. CONCLUSIONS 

The experiment reported here provides a highly reproducible violation of quantum me- 
chanical duality using an apparatus configured with two independent HeNe lasers and readily 
available components. The determination of that duality violation resolves to easily mea- 
sured variations in continuous wave laser beam power, variations predicted by a locally real 
representation of quantum mechanics but excluded by the probabilistic interpretation of 
quantum mechanics. 
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